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ANALYTICAL EXPRESSIONS FOR POSITION ERROR IN 
TRIANGULATION SOLUTION OF POINT IN SPACE 
FOR SEVERAL STATION CONFIGURATIONS 

By Sheila Ann T. Long 
Langley Research Center 

SUMMARY 

Analytical expressions are derived to first order for the rms position error in the 
triangulation solution of a point object in space for several ideal observation- station con- 
figurations. These expressions provide insights into the nature of the dependence of the 
rms position error on certain of the experimental parameters involved. The station 
geometries examined are: (1) the configuration of two arbitrarily located stations; (2) the 
symmetrical circular configuration of two or more stations with equal elevation angles; 
and (3) the circular configuration of more than two stations with equal elevation angles, 
when one of the stations is permitted to drift around the circle from its position of sym- 
metry. The expressions for the rms position error are expressed as functions of the 
rms line-of-sight errors, the total number of stations of interest, and the elevation angles. 

INTRODUCTION 

Obtaining the errors in the triangulation solution of a point object or an elongated 
object in space using data from three or more arbitrarily located observation stations 
is a complicated problem, and numerical solutions are usually sought. Reference 1 gives 
formulas for the errors in the geocentric position, as calculated using the simultaneity 
circle, of a satellite. Reference 2 presents formulas for adjusting all measurements of 
a space triangulation to determine the coordinates of the stations on the basis of two or 
more well-known stations. In reference 3 formulas are derived for the rms errors in 
location, orientation, and shape in the triangulation solution of an elongated object. By 
programing the formulas in references 1, 2, and 3 for a high-speed computer and using 
data for specific object- station relationships, numerical results can be computed. Ref- 
erences 4, 5, and 6 give numerical results for the following respective satellite triangu- 
lation nets: five specific stations in the United States; nine specific stations covering the 
Australian continent; and thirty- six specific stations covering the earth. 

Numerical studies, however, possess the limitation that they do not provide com- 
plete insights into the nature of the dependence of the errors on the experimental parame- 



ters involved. Several analytical studies which do provide some of these insights can 
be found in the literature. In reference 7 analytical expressions for the position error 
in the solution of a point in space are presented for a system employing angle-only infor- 
mation for two stations and a system employing range -only information for three stations. 
For this latter system the three-station configuration is separated into three two-station 
configurations and then the separate results are combined to give the total result. Ref- 
erence 8 gives analytical expressions for the position error in the range trilateration 
solution of a satellite for two ideal station configurations. The first configuration has 
three stations arranged such that the angles 3,1,2 and 1,2,3 are equal. The second con- 
figuration has four stations arranged such that the angles 3,1,2; 1,2,3; 2,4,3; and 2,3,4 are 
all equal. For both configurations the satellite is located at the zenith as referenced to 
the centroid of the triangle 1,2,3. 

The purpose of the present paper is to derive analytical expressions for the rms 
position error in the triangulation solution of a point in space for several additional ideal 
observation- station configurations. These expressions are to provide additional insights 
into the nature of the dependence of the rms position error on certain of the parameters 
involved. 

The first observation- station configuration examined is the configuration of two 
arbitrarily located stations. Since, in general, the lines of sight from two stations do 
not intersect, a most probable point in space must be chosen. Two choices for the most 
probable point - the midpoint of the shortest line between the two lines of sight and the 
point, on the aforementioned line, which subtends equal residual angles at the two 
stations - are discussed and compared. The situation in which the data from one of 
the stations is degraded relative to that from the other is considered. For the first 
choice of the most probable point, an expression for the optimum relative weighting 
factor for the data from the first station, which minimizes the rms position error, is 
derived. 

The second observation- station configuration examined is the symmetrical circu- 
lar configuration of two or more stations with equal elevation angles. The third configu- 
ration examined is identical to the second configuration, except that there are more than 
two stations and one of them is permitted to drift around the circle from its position of 
symmetry. 

For each observation- station configuration considered, an analytical expression for 
the rms position error in the triangulation solution of a point in space — as a function of 
the rms line-of-sight errors, the total number of stations of interest, and the elevation 
angles — is derived to first order. Also for each configuration considered, the optimum 
elevation angle, which minimizes the rms position error, is found. Considerations are 
given as to whether the symmetrical arrangement of the circular configuration of more 
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than two stations with equal elevation angles is the arrangement of this configuration 
which yields the minimum rms position error. 
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SYMBOLS 

coefficients depending on n, y, and 9 
actual point object in space 

distance from each station in circular configuration to point A g 
distances from stations Sj and S 2 , respectively, to point A g 
distance between stations in two- station configuration 

line-of-sight errors in plane of 9 from pth and qth stations, respectively 

line-of-sight errors in plane of 9^ and f rom stations Sj and S 2 , 
respectively 

line-of-sight errors out of plane of 9 from pth and qth stations, respectively 

line-of-sight errors out of plane of 9^ and 9 2 from stations and S 2 , 
respectively 

position error function 

line-of-sight unit vector, in direction of line of sight from pth station 
coefficients depending on A,B,C,D,E,F, and jd 
altitude of point A g 

actual unit vector, in direction from pth station to point A g 
point defined on figure 1 

total number of stations in circular configuration 
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most probable point in space 

indices for sequential labeling of stations in circular configuration 

radius of circular configuration of stations 

most probable vector, from pth station to point P s 

first and second stations, respectively, in configuration of two arbitrarily 
located stations 

symmetrical determinant of coefficients A,B,C,D,E,F 
sin2e 

sin 0^ and sin 0 2 , respectively 

rectangular coordinate system with origin at center of circular configuration 
of stations, x-axis toward original position of Oth station, z-axis toward 
point A s , and y-axis to form right-hand orthogonal triad 

unit vectors in directions of increasing x, y, and z, respectively 

angles which pth and qth stations, respectively, make with x-axis 

variable angle which drifting station makes with x-axis 

difference vector, between unit vector ip and unit vector Cp 

Krone eke r delta 

x-, y-, and z- components, respectively, of vector 6 r 

line-of-sight error vector, shortest vector from line of sight from pth 
station to point P g 

position-error vector, from point A g to point P g 

rms position error in triangulation solution of point in space, in general 
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e* rms position error in triangulation solution of point in space - for two arbi- 

trarily located stations with point P g taken as point, on shortest line 
between lines of sight, which subtends equal residual angles at stations 

£,?),£ coefficients depending on a, d0 , d0 , a , and Q 

r r P 

8 elevation angle, in general 

elevation angles from stations and Sg, respectively 

a rms line-of-sight error, in general 

rms line-of-sight errors from stations Sj and S2, respectively, when data 
from one station is degraded relative to that from the other 

w relative weighting factor for data from station 

A bar over an expression is used to denote the mean value. 

ANALYTICAL FORMULATION 

Position Error in Solution of Point for 
Two Arbitrarily Located Stations 

The two arbitrarily located observation stations are denoted by and S2, 
respectively, and are shown in figure 1 . The two stations are separated by the distance 
b. The angles from the stations Sj and S2 to the actual point object A g in space 
are denoted by and $2, respectively; these angles will be herein called elevation 
angles. The distances of the point A g from the stations and S2 are denoted by 
a^ and a^, respectively. The distance from the base line to the point A s is denoted 
by h; this distance will be herein called the altitude. In practice, errors, herein called 
line-of-sight errors, will exist in the measured values of the elevation angles 8 ^ and 
6?2- The line-of-sight errors in the plane of the angles 6 ^ and 8 2 are denoted by d 0 ^ 
and d02> respectively; the line-of-sight errors out of the plane of the angles and 62 
are denoted by d0^ and d0g, respectively. 

The position error function e in the triangulation solution of a point in space for 
two arbitrarily located observation stations is given by 

de = -|f- d0 1 + d0 + d 0„ 

d 8 1 1 90 2 1 9 0 1 1 90 2 2 
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The rms position error e in the solution of a point for two arbitrarily located stations is 
then given by 



d0, +|^-d0, + -^-d0 1 + 
1 90 2 1 d <P 1 1 



( 1 ) 


where the bar denotes the mean value. It can be assumed that the line-of -sight errors 
d0p d$ 2 , d0-p and d0 2 are uncorrelated; hence, equation (1) becomes 



Since each of the errors da^, d$ 2 , d0j, and d 0 2 is perpendicular to its line of sight, 
it can be assumed that 


d flj 2 = d6 l 2 2 = d^ 2 = d0 2 2 = a2 (3) 

where a is the rms line-of-sight error. Therefore, equation (2) becomes 


e2 = 0-2 


Se 

dG< 


3e 

30, 


3e 

d( P i 


3e 

30 2 


(4) 


From figure 1, using the law of sines yields 

a l _ b 

sin 0 2 sin^ + 9<^ 

Hence, the distance a^ of the actual point A g in space from the observation station 
b sin $2 

3^ = 


is 


1 sin^ + 0 2 ^ 


(5) 


Similarly, the distance a 2 of the point A s from the station S 2 is 
b sin 6 - l 


a 2 = ~ 


sin 


(Si + e 2 ) 


(6) 


The line-of-sight errors d0p d0 2 , d0^, and d0 2 a.re all independent and, hence, 
their respective effects on the rms position error in the triangulation solution of a point 
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in space can be examined separately. The position error resulting purely from the error 
is examined first. From figure 1 it is observed that the error d<9^ manifests itself 
in a change in the distance ^ > while the elevation angle 8^ * s unaffected. Hence, 

Be _ 9a 2 _ k sin 6 ^ 

9 ^i sin^gj + 

Similarly, 

_0e_ = ^4 = b sin h (g) 

80 2 9f? 2 sin 2 ^ + ^ 2 ) 


Next, the position error resulting purely from the line- of- sight error d0^ is 
examined. The error d0^ is out of the plane of the elevation angle 0^. This pre- 
vents the lines of sight from the observation stations from intersecting. For this situa- 
tion a most probable point P g in space must be chosen. Two choices are treated in the 
present paper. 

Most probable point taken as midpoint of shortest line between lines of sight .- One 
choice for the most probable point P s in space is the midpoint of the shortest line 
between the two lines of sight. This choice is used in reference 9. From figure 1 it is 
seen that the line segment A g P g is 


A S P 


s 


A S N 

2 


Since dcp^ is small, 

A s P s = 2 a l d ^l 


( 9 ) 


Hence, 


Be 



Using equation (5) yields 


3e 


B ( p 1 


b sin 8 2 
2 sin^ + # 2 ) 


Similarly, 


( 10 ) 


7 



8e 


(id 


b sin 

9 $2 2 sin^0j + ^2) 

Substituting equations (7), (8), (10), and (11) into equation (4) produces 

^ (sin 2 ^ + sin20 2 jfl + ± Sin2(0 1 + 0 2 )| • (12) 

e 2 - a 2b 2 1 ■—= 1 

sin**^ + 0 2 ) 


Figure 1 shows that the distance b between the two observation stations is 

b = cos 0^ + a 2 cos # 2 (13) 

where the distances ^ and a 2 of the actual point A s in space from the stations Sj. 
and S 2 , respectively, are 


a 1 = — — 
1 sin 0^ 


(14) 


a 9 = — 0- (15) 

z sin 0 2 

where h is the altitude of the point A s . Substituting equations (14) and (15) into equa- 
tion (13) leads to 


b = 


h sin(0j_ + 0 2 ) 
sin Q 1 sin # 2 


Substituting equation (16) into equation (12) results in 


e 2 = a 2 h 2 


^5 in 2 #! + sin 2 # 2 j 

1 + | sin2^ x + # 2 ) 

sin 2 #! sin 2 #. 

2 + e 2 ) 


(16/ 


(17) 


Therefore, the rms position error e in the triangulation solution of a point in 
space for two arbitrarily located observation stations, with the most probable point in 
space taken as the midpoint of the shortest line between the two lines of sight, is 
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e - ha < 


r~ 

(sin 2 ©^ 


+ sin 2 0 2 j|l + sin 2 ^ 0 . 


1/2 


H) 


sin 20 ^ sin 2 0 2 sin 2 ( 0 i + 02 ) 


(18) 


It is noted that the rms position error is directly proportional to the rms line-of- sight 
error. It is also noted that the rms position error is infinite for 0^ + 0 2 = 180°, which 
occurs when the two lines of sight are parallel. 

Table' I is a table of values of e/ha as a function of the elevation angles 0-^ and 
02* For convenience, the only values of e/ha listed are those for the angles 0^ and 
02 at intervals of 10°. From table I it is seen that the minimum value of e/ha occurs 
for 6-y equal to 02 In the neighborhood of 55°. 

Then, the quantity e/ha was evaluated for values of 0^ and 02 from 45° to 65° 
at increments of 0.001°. Table II is a table of values of e/ha as a function of 0^ and 
02 at increments of 0.001°. For convenience, the only values of e/ha listed are those 
for 0^ and 0 2 from 56.096° to 56.100°. From table II it is observed that the minimum 
value of e/ha occurs for 0^ = 02 = 56.098°. Therefore, the minimum rms position 
error in the triangulation solution of a point in space for two arbitrarily located obser- 
vation stations, with the most probable point in space taken as the midpoint of the shortest 
line between the two lines of sight, occurs when the two elevation angles are such that 
0 X = 0 2 = 56.098°. 

Setting 0j = 0 2 = 0 in equation (18) establishes the following: 


oh /l + sin 2 0 - sin^ ^ , 1Q v 

“ M sin 4 0 - sinS ) ‘ ’ 

Therefore, equation (19) is the equation for the rms position error in the triangulation 
solution of a point in space for two arbitrarily located observation stations, with the most 
probable point in space taken as the midpoint of the shortest line between the two lines of 
sight, when the two elevation angles are set equal. 

Most probable point taken as point, on shortest line between lines of sight, which 
subtends equal residual angles at stations .- A second choice for the most probable point 
P g in space is the point, on the shortest line between the two lines of sight, which sub- 
tends equal residual angles at the two observation stations. This choice is used in ref- 
erence 8. From figure 1 it is seen that 
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Or, 


A S P S __ A g N - AgP s 
a 2 " a l 


As p s 


a 2 

( a l + a 2) 


A S N 


Since d <p^ is small, 


A P = 


a^agdc^ ^ 

( a l + a 2) 


( 20 ) 


Hence, 


ae _ a i a 2 
90! ” ( a l +a 2 ) 


Using equations (5) and (6) results in 


8e 
9 4>i 


b sin 0^ sin 6 2 
(sin 61 + sin 02) sin (#i + #2) 


Because of symmetry 


( 21 ) 


ae 

902 


b sin 6 1 sin 0 2 


(sin + sin 0 2 jsin(0j + 0 2 ) 


( 22 ) 


Changing the notation from e to e' in equation (4) and then substituting equa- 
tions (7), (8), (21), and (22) into the equation for (e') 2 result in 


(e’)2 - <72)32 


2 sin2^ 1 sin 2 0 2 sin 2 (0-) + 0 2 ) 
Sln8 l + sin e 2 + r~ — —~To 


sin' 


(sin 9\ + sin 6)2)' 
+ 6 2) 


(23) 


Therefore, equation (23) is the equation which gives the rms position error in the triangu- 
lation solution of a point in space for two arbitrarily located observation stations - with 
the most probable point in space taken as the point, on the shortest line between the two 
lines of sight, which subtends equal residual angles at the stations. 
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Comparison of position errors for two choices of most probable point .- Taking the 
ratio of equation (23) to equation (12) renders 


sin 2 ^! + 0 2 ) 

ijsin 2 ^ + sin 2 0 2 )( sin + sin ^) 2 “ 2 sin20 l sin 2 0 2 

(sin 0^ + sin 0 2 ) 2 ^sin 2 0^ + sin 2 £ 2 j 

jl + | sin 2 (0 1 + 0 2 ) 


Hence, 

~ = ^ for ijsin 2 ^ + sin 2 0 2 )(sin + sin 0 2 )^ - 2 sin 2 ^ sin 2 0 2 = 0 


(24) 


(25) 


For convenience, the following change of variables can be made: Xj = sin 
X 2 = sin 02- Since the elevation angles and 0 2 are positive, then the variables 
Xj^ and X 2 are also positive. Making this change of variables in expression (25) 
furnishes 


jSl for i(x x 2 + X 2 2)(X 1 + x 2 ) 2 - 2 Xj2x 2 2 g 0 

If Xj = X 2 , then jx^ 2 + X 2 ^) = 2 XjX 2 and the following result can be established: 

7=1, since i(Xj 2 + X-^Xj + X 2 ) 2 - 2 Xl 2x 2 2 = 0 

If X^ * X 2? then ^X^ 2 + X 2 2 ) > 2X^X 2 and the following result can be proved: 

e i < 1, since |(X! 2 + X 2 2 )(X 1 + X 2 ) 2 - 2 Xl 2x 2 2 > 0 
Consequently, 

e r = e /$ 1 = 9 2 \ (26) 

e' < e ( 0 1* 0 2) ( 2 ?) 

Hence, when the two elevation angles are equal, the rms position error e’ is equal 
to the rms position error e. However, when the two elevation angles are not equal, then 
€' is less than e. Therefore, the choice of the most probable point in space as the point, 
on the shortest line between the two lines of sight, which subtends equal residual angles 
at the two observation stations, yields a smaller value of the rms position error than does 
the choice of the most probable point as the midpoint of the aforementioned line, except 
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in the special situation of equal elevation angles and then the two choices yield the same 
value . 

Optimum relative weighting factor for data from first of two arbitrarily located 
stations.- Earlier in this analysis f or tw o arbit ra rily lo cated obser vation stations, the 
mean- square line -of -sight errors d^ 2 , d0 2 2 , dc^ 2 , and d<£ 2 2 were assumed to be 
all equal. The situation in which the mean- square line-of- sight errors from the station 
S-^ differ from the mean- square line -of- sight errors from the station S 2 are now con- 


Hence, 


dO^ = d$^ 2 - Oj 2 

(28) 

d0 2 2 = d< # , 2 2 = °2 2 

(29) 


where Oj and a 2 are the rms line-of- sight errors from the stations and S 2 , 
respectively. Substituting equations (28) and (29) where appropriate into equation (2) 
gives 


e 2 = 



(30) 


Since the rms line-of-sight errors from the two observation stations are different, 
the data from the two stations must be weighted differently. The relative weighting fac- 
tor for the data from the station can be denoted by t o. The purpose of this section 
is to determine the optimum value of the relative weighting factor w, which minimizes 
the rms position error in the triangulation solution of a point in space. 

For a line-of-sight error d$j at the observation station S^, shown in figure 1, 
the relative weighting factor cu is 

A s P s 

" = A s N - A s p s 


Or, 


A S P 



A S N 


Since d<p is small, 
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Hence, 


A s p s = 


(31) 


3e = / o> \ 3 
90^ \1 + ay 1 


Using equation (5) renders 


a D , b sin 6 o 

<*e / to \ 

9 \1 + o>/sin^ + Q^j 


For an error d<£ 2 at the station S 2 , the relative weighting factor a> is 


S N A S P S 
A S P 


(32) 


Or, 


A P = 

As s u + ur s 


IA.N 


Since d<^ 2 Is small, 

A s p s = (iT^ a 2 d< #>2 

i 

Hence, < 

Using equation (6) gives 
9e 


9 cf> t 


1 \ b sin 9 1 


1 + (j 0 / 


sin^0j + 0 2 ) 


Substituting equations (7), (8), (32), and (34) into equation (30) furnishes 
sin^e^ + ^<t>2 a ^2 s i n 2g^ + a 2 sln^gjJ 


e2 = b 2 


+ 0 2 ) 


(1 + o>) Sin2^0j + 


(33) 


(34) 


(35) 
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In order to find the optimum relative weighting factor, the partial derivative of equa- 
tion (35) with respect to o> is taken and then this resulting expression is set equal to 
zero. That is, 


00) 


(1 + cu)^ sin 2 ^ + 02)(2 cl>ctj 2 gin 2 ^) 


(1 + our sin 4 (0, + e 2 ) 


(a; 2 sin 2 0 2 + a 2 2 sin20 l 

^(2 + 2cu)sin2 

( 0 i + h) 

(1 + a.) 4 sln 4 ( 

B 1 + B z) 

_ 


= 0 


(36) 


Hence, 

cu(l + cu)ct^ 2 sin2f? 2 = ou^cr^ sin^Sg + 
Consequently, 

a 0 2 sin 2 0i 

w ~ ~^2 2 

<7^ sin^f?2 


(37) 


Therefore, for the choice of the most probable point in space as the midpoint of 
the shortest line between the two lines of sight, equation (37) is the equation for the 
optimum relative weighting factor, which minimizes the rms position error in the tri- 
angulation solution of a point in space, for the data from the observation station S^ 
when the data from one of the stations is degraded relative to that from the other. It is 
noted from equation (37) that the optimum relative weighting factor for the data from the 
station S-^ is directly proportional to the square of the ratio of the rms line -of- sight 
errors from the stations S 2 and S-^. 


Position Error in Solution of Point for Symmetrical 
Circular Configuration of n g 2 Stations 
With Equal Elevation Angles 

The configuration of n, where n t. 2, observation stations is now examined. The 
n stations are equally spaced around the circumference of a circle, the center of which 
is located at the foot of the perpendicular drawn from the actual point A s in space, as 
shown in figure 2. For this configuration an approach different from that of the preceding 
section is used. 
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The actual unit vector i p is in the direction from the pth observation station to the 
actual point A s in space. The distance between each of the n stations and the point 
A s is denoted by a. The most probable vector r p is the vector from the pth station 
to the most probable point P g in space. The point P g is displaced from the point A s 
by the position error vector F r . The line of sight from the pth station is in the direction 
of the line-of-sight unit vector e p . The unit vector e p is displaced from the unit vec- 
tor ip by the difference vector A p . Lastly, the line-of-sight error vector, the shortest 
vector from the line of sight from the pth station to the point Pg, is denoted by 6 p . 

From figure 2 it is observed that the line-of-sight error vector 5 p , the most 
probable vector r p , and the line-of-sight unit vector e p , respectively, from the pth 
observation station are 



r p = ai p + 6 r 


(39) 


- *p + A P 


(40) 


Substituting equations (39) and (40) into equation (38) renders 

= aL + fiL - ( aL + 6 * ( L + A t 


'p - -r u r - (ai p + 6 r ) * (i p + ^ p )(ip + *p) 


Keeping only terms to first order and noting that i p ■ A p = 0 to first order produce 


5 p = 6 r 


- (‘P ' «r) 


^ ' aA p 


Hence, again noticing that ip * A p = 0 to first order leads to 

(«p ) 2 = («r ) 2 - (*p ' 6 r ) 2 + a2 (^p ) 2 - 2a(6 r • S p ) 

Therefore, for the n stations equation (42) becomes 

n-1 n-1 n-1 _ n-1 


n-1 


(41) 


(42) 


X (®p) 2 = I (®r) 2 ' I ( f P • ®r) 2 + a2 I (*p) 2 - £ (*r ' *p) <«> 

p=0 p=0 p=0 p=0 p=0 


where n = 2. 
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The elevation angle from each of the n observation stations is the same, since the 
n stations are located on the circumference of a circle with the actual point A s in 
space being located on the perpendicular whose foot is at the center of the circle. The 
elevation angle is denoted by 9 and is shown in figure 2. The angles d0 p and d$ p 
are the line -of -sight errors in and out of the plane of the angle 9, respectively, from the 
pth station. The angle a p is the angle between the x-axis and the pth station, and 
rv = 2-rm/n whpro n = 0. 1 , 2, . , n - 1 , as the n stations are equally spaced around 

w '- p — ■■£'(' "7 ” a. 7 7 7 ? j *■ 

the circumference of the circle. The rectangular coordinate system of figure 2 has its 
origin at the center of the circular configuration of stations, x-axis toward the Oth station, 
z-axis toward the point A s , and y-axis such as to form a right-hand orthogonal triad. 

In order to determine the rms position error in the triangulation solution of a point in 
space as a function of the angle 8 for this configuration, equation (43) must first be 
written as a function of the angle 0, the errors d0 p and d^ p , the angle o p , and the 
x-, y- ? and z-components of the position error vector 6 p . 

The actual unit vector ip from the pth observation station corresponds to the sit- 
uation in which no line -of- sight errors exist in the line of sight from the pth station (l.e., 

\ A \ 

when d0 p = d0 p = 01. The unit vector i p from the pth station is 


ip = -x cos 9 cos a p - y cos 8 sin a n + z sin 8 


(44) 


where x, y, and z are the unit vectors in the directions of increasing x, y, and z, 
respectively. From figure 2 it is observed that the difference vector A p from the pth 
station is 


= e p - i p (45) 

where e p is the line-of- sight unit vector from the pth station. Hence, 

A p = x^d0 p sin 6 cos ot p - d<£ p sin a p j + y^d0 p sin 9 sin Q! p + d<£> p cos a p j + z d0 p cos 6 

(46) 

The position error vector ff r can be expressed as 

5 r - x6 x + y6y + z6 z (47) 

where 5 X , 6y, and 6 Z are its x-, y-, and z-components, respectively. 

Using equations (44), (46), and (47) establishes the desired functional dependence of 
each term in equation (43). Hence, 
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(48) 


n-1 n-1 


I (^r) 2 = £ ( 6 x 2 + 6 y 2 + 6 z 2 ) = n ( 6 x 2 + 6 y 2 + 5 z 2 ) 

P-0 p=0 


n-1 2 2 

Z (‘p • ■ I (a cos 6 cos on p - 6y cos 9 sin a p + 6 Z sin e'j 

p=0 p=0 

n-1 n-1 

= 6 X ^ cos2fl ^ cos2ci?p + 6y2 cos2fl sin2o! p + n6 z ^ sin2# 
p=0 P=0 

n-1 n-1 

- 26y5 z cos 9 sin 9 ^ sin otp - 26 X 6 Z cos 6 sin 0 ^ cos a p 

p=0 p=0 

n-1 

+ cos2# ^ sin^QTpj (49) 

P=0 


n-1 n-1 

^ (^p ) 2 = a2 2, |( d0 p sin 6 cos a p " d ^p sin a p ) 2 

p=0 p=0 L 

+ ^d0p sin 6 sin a p + d<£ p cos + (d0 p cos (fj 


(50) 


n-1 n-1 

2a ^ ■ A p ) = 2a6 x ^ (dd p sin 6 cos a p - d0 p sin a 

p=0 p-0 

n-1 

+ 2a 6 y ^ ^d0p sin 0 sin a p + d^ p cos 

p-0 

n-1 

+ 2a 6 Z cos 0 d0 p (51) 

P=0 


Substituting equations (Al), (A2), (A3), (A4), and (A5), which are derived in the appendix, 
into equation (49) produces 
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n -1 _ 6 2 cos 2 e + 0 : 6 2 cos 20 + n g 2 S in 2 g 

I(f p -a r ) 2 = 2 

p = 0 l 25 x 2 cos20 + 26 z 2 sin 2 0 


(n > 2 ) 


(n = 2 ) 


( 52 ) 


The substitution of equation (52) into equation (43) yields two different results, one 
for n > 2 and one for n = 2. Hence, these two situations must be examined separately. 

Symmetrical circular configuration of n > 2 stations with equal elevation angles 
Substituting equation (52) for n > 2 along with equations (48), (50), and (51) into equa- 

n - 1 2 

tion ( 43 ) leads to an expression for ^ for n > 2 as a function of the elevation 

p =0 

angle 6 , the line-of-sight errors d 0 p and d<£ p , the angle a p , and the x-, y-, and 
z- components of the position error vector 6 r - That is, 


n -1 


1 (®p) = "(' 

p=0 


5 2 + 5 2 + 6 2 

x y z 




fn K 2 , n K 2 4 . n * 2 

1 — W6“P T “ T - UU„ Lli.ll V/ ) 


+ a< 


n -1 

t | 

p=0 ^ 


2 


d 0 sin 6 cos a n - sin 

F 


7 


p “V'p 

2 


+ ^d 0 p sin 9 sin a p + d 0 p cos a^y + ^d£? p cos By J 
n -1 

- 2 a 6 x ^ ^da p sin 9 cos a p - d$ p sin 


P=0 

n -1 


- 2aSy ^ ^d0 p sin 9 sin a p + d<£ p cos a p j 

p =0 

n -1 

^ d 9 n (n > 2 ) 

p =0 


- 2 a 6 z cos 6 / dfl p 


( 53 ) 


n -1 


In this form, can be minimized with respect to the x-, y-, and 

p=0 

z- components of the vector ff r . 
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To obtain the minimum with respect to 6^ the partial derivative of equation (53) 
with respect to 6 X is taken and then this resulting expression is set equal to zero. That 
is, 


Or, 



n-1 

2n6 x - n6 x cos 2 # - 2a Y /d# p sin # cos a p - d^ p sin a p ) = 0 

p=0 7 


n-1 

(l + sin 2 #)n6 x = 2a ^ (d # p sin 9 cos a p - d^ p sin a^j 

p=o V 


(54) 


Thus, 


n-1 n-1 


(l + sin 2 #) n 2 6 x 2 = 4a 2 ^ ^ /d# p sin 0 cos a p - d <£ p sin a p Vd# q sin 9 cos a q - d<£ q sin a q \ 

p=0 q=0 ' 


(55) 


where the bars denote the mean values. 


It can be assumed that all of the line-of- sight errors d# p and d^ p are uncorre- 
lated. Also, since the errors d# p and d<£ p are all perpendicular to their respective 
lines of sight, it can be assumed that 


d9 p d0 q = ^Spq 


(56) 


a ^pq 

(57) 

d# p d(^)q — d#qd^ p = 0 

(58) 


where 6 pq is the Krone cker delta and a is the rms line -of -sight error. 

Substituting equations (56) to (58) where appropriate into equation (55) results in 


n-1 


n-1 


(l + sin 2 #)^n 2 6 x 2 = 4a 2 o 2 |sin 2 # ^ cos 2 a p + ^ sin 2 apj 

p=0 p=0 


V 


/ 


Substituting equations (A4) and (A 5) for n > 2 into equation (59) gives 

\2 


(l + sin 2 #) n 2 6 x 2 = 4a 2 cr 2 ^ sin 2 # + 


(59) 
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Therefore, 


5~2 = 2a 2 cr 2 

n(l + sin 2 #) 


( 60 ) 


Minimizing equation (53) with respect to 6y and 5 Z , respectively, results in the 
following two en nations: 


n-1 


(l + sin 2 #)n6 y = 2a ^ /d# p sin 8 sin a p + d</) p cos a p \ 
p=o v 


(61) 


n-1 


P=0 


n£> z cos 9 = a J d# p 


(62) 


Using equations (61) and (62), respectively, and procedures analogous to those used for 
obtaining equation (60) for & x 2 yields the following equations for 6 y 2 and 6 Z 2 : 


6, 2 = 


2a 2 a 2 


n(l + sin 2 #) 


(63) 


i2rr2 


2 _ ^ g 
°z “ *> 

n cos 0 


(64) 


The rms position error e in the triangulation solution of a point in space for the 
circular configuration of n observation stations with equal elevation angles is given by 


e 2 = 6 X 2 + 6 y 2 + 6 Z 2 


(65) 


Substituting equations (60), (63), and (64) into equation (65) furnishes 


5 _ a 2 o 2 /s - 3 sin 2 #' 
6 = " l- S in 4 J 


( 66 ) 


From figure 2 it is seen that the distance a from each of the n stations to the actual 
point A g in space is 


a = 


sin 8 


(67) 
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Substituting equation (67) into equation ( 66 ) renders 


e 2 - 


h 2 0 r 2 / 5 _ 3 sin 2 0 


n 


(Sin^0 - sin®0y 


( 68 ) 


Therefore, the rms position error e in the triangulation solution of a point in 
space for the symmetrical circular configuration of n > 2 observation stations with 
equal elevation angles is 


e 


ho/ 5-3 sin 3 0 
^\sin 2 0 - sin^ 0 , 


1/2 


(n > 2 ) (69) 


The rms position error from equation (69) is directly proportional to the rms line- of - 
sight error and inversely proportional to the square root of the total number of stations 
of interest. For elevation angles of 90° the rms position error is infinite. For an ele- 
vation angle of 90° the lines of sight from the different stations could not intersect. Fig- 
ure 3 is a plot of e/ha for n > 2 as a function of the elevation angle 0 , for n = 3, 4, 

5, 6 , 7, and 8 . 

The optimum elevation angle, which minimizes the rms position error in the triangu- 
lation solution of a point in space, for the symmetrical circular configuration of n > 2 
observation stations with equal elevation angles can be determined by minimizing equa- 
tion ( 68 ) with respect to the elevation angle 0 . That is, 


9e 2 _ 2 h 2 o 2 sin 0 cos 0 

-6 sin ®0 + 15 sin^0 - 5 

90 n 


^sin 2 0 - sin® 0 , 

2 


Hence, 

-6 sin ®0 + 15 sin ^0 -5 = 0 

Making the change of variables X = sin 20 in equation (71) gives 
X 3 - | X 2 + | = 0 

The three roots of this equation are 

X = 0.675905; -0.524875; 2.348970 


(70) 


(71) 


(72) 
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Since X = sin^fl, the root corresponding to the physical solution must be positive and 
must be between zero and unity. Hence, the root corresponding to the physical solution 
is X = 0.675905. Therefore, 


9 = sin-iX 1 / 2 = 55.300° (n > 2) (73) 

Thus, the optimum elevation angle, which minimizes the rms position error in the solution 
of a point, for the symmetrical circular configuration of n > 2 stations with equal ele- 
vation angles is 55.300°. This result is also seen from figure 3. In addition, for 
9 = 55.300° ± 5° (see fig. 3), the difference in the rms position error for n = 3 is 
2.5 percent, while for n = 8 it is 0.9 percent. 

The optimum elevation angle (55.300°) is within a fraction of a degree of the angle 
(54.667°) which the sloping edge of a tetrahedron makes with its base. Therefore, the 
optimum configuration of three observation stations for triangulating on a point in space 
is the configuration of the three stations plus the point which closely approximates a 
tetrahedron. 

Given the altitude of the point in space, the value of the optimum elevation angle for 
n > 2 can be used to compute the value of the optimum radius of the symmetrical circular 
« ri «f4 £ y 1irr iHoTi of n > 2 observation stations with equal elevation angles. The radius R 
of the circular configuration (see fig. 2) is 

R = h cot S (^4) 

Hence, substituting 9 = 55.300° into equation (74) produces 

R = (0.69243)h (n > 2) (75) 

Therefore, equation (75) furnishes the value of the optimum radius of the symmetrical 
circular configuration of n > 2 stations with equal elevation angles, for a given altitude 
of the point. 

Symmetrical circular configuration of n — 2 stations with equal elevation angles . - 
Noticing that a p -up for n = 2 and substituting equations (48), (50), (51), and (52) for 

1 2 

n * 2 into equation (43) establish an expression for ^ ^p) as a function of the ele- 

p=0 

vation angle 0, the line-of- sight errors d0 p and d^ p , the angle a p * ?rp, and the x-, 
y-, and z-components of the position error vector 6 r . That is, 
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^ ( 6 p ) 2 = 2(a x 2 + 5 y 2 + 6 Z 2 ) - (26 x 2 cos 2 0 + 25 z 2 sin 2 s) 
p=o 1 

+ a 2 ^ j|d0p sin 9 cos(ttp ) - d<£ p sin(7ip)j 2 


p=0 


+ ^de p sin 0 sin(Trp) + d<£ p cos(np)j 2 + ^d0 p cos 0j 2 
1 

- 2a6 x ^ jd0 p sin 9 cos(ttp) - d$ p sin(7rp)j 


P-0 

1 


- 2a 6y p0p sin 9 sin(?rp) + d0 p cos(ttp)] 

p=0 L 


- 2a6 z cos 9 / d0 p 


P=0 


(n= 2) 


(76) 


Using procedures analogous to those used in arriving at equation (68) from equation (53) 
yields 


2 _ h 2 o~ 2 (\ + sin 2 g - sin^ 9 ^ 


e“ = 


sin^f? - sin®0 


(77) 


Therefore, the rms position error e in the triangulation solution of a point in 
space for the symmetrical circular configuration of n = 2 observation stations with 
equal elevation angles is 

vl/2 

£ = hg/l + sin^slnM (n = 2 ) (78) 

sin 4 S - sin® 8 j 

This equation is identical to equation (19), which gives the rms position error in the solu- 
tion of a point for two arbitrarily located stations, with the most probable point taken as 
the midpoint of the shortest line between the two lines of sight, when the two elevation 
angles are set equal. 
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From equation (78) it is noticed that the rms position error in the triangulation 
solution of a point in space for the symmetrical circular configuration of n = 2, as for 
n > 2, observation stations with equal elevation angles is directly proportional to the rms 
line -of -sight error and infinite for elevation angles of 90°. Figure 4 is a plot of e/ha 
for n = 2 as a function of the elevation angle 9. 

The optimum elevation angle, which minimizes the rms position error in the trian- 
guiation solution of a point in space, for the symmetrical circular configuration of n = 2 
observation stations with equal elevation angles can be determined by minimizing equa- 
tion (77) with respect to the elevation angle 9. That is, 


a£ = h 2 CT 2 sin 3 0 cos e f( sln2e - f^Xj - 2 sin2g ) - (j + sin2g - Sln4s )( 2 - 3 Sln2 s)| = 


a 9 


Hence, 


(sin*0 - sin**#) 


(79) 


(sin 3 ^ - sin4fl)(l - 2 sin2fl) - (l + sin2# - sin4#^2 - 3 sin2^ = 0 (80) 

Making the change of variables X = sin2# in equation (80) renders 

X 3 - 2X 2 - 2X + 2 = 0 (81) 


The root corresponding to the physical solution is X = 0.688892. Therefore, 


6 a sin" ix 1 / 2 = 56.098° (n = 2) (82) 

Thus, the optimum elevation angle, which minimizes the rms position error in the solu- 
tion of a point, for the symmetrical circular configuration of n - 2 stations with equal 
elevation angles is 56.098°. This result is also seen from figure 4. In addition, for 
9 = 56.098° ± 5° (see fig. 4), the difference in the rms position error for n = 2 is 
1.5 percent. 

It is recalled that the minimum rms position error in the triangulation solution of 
a point in space for two arbitrarily located observation stations, with the most probable 
point taken as the midpoint of the shortest line between the two lines of sight, occurs 
when the two elevation angles 6 j and 0 2 are such that ^ = #2 = 56 * 09 8°. Hence, a 
necessary and sufficient condition for a minimum rms position error in the solution of a 
point for two arbitrarily located stations is that the two elevation angles be equal, and 
equal to 56.098°. 

Substituting 9 = 56.098° into equation (74) renders 

R = (0.67200)h (n = 2) (83) 
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Therefore, equation (83) furnishes the value of the optimum radius of the symmetrical 
circular configuration of n = 2 observation stations with equal elevation angles, given 
the altitude of the point in space. Since the distance b between the two stations is twice 
the radius R of the circular configuration, then 

b = 2R = (1.3440)h (84) 

Therefore, equation (84) produces the value of the optimum distance between two stations 
with equal elevation angles, given the altitude of the point. 

Variation in Position Error in Solution of Point for Circular Configuration 
of n > 2 Stations with Equal Elevation Angles, When One 
Station is Permitted to Drift Around the Circle 

For this part of the analysis, the Oth observation station ^i.e., the one, in the previ- 
ous part, that was located on the x-axis and, hence, for which = «q * 0) is permitted 
to drift around the circle. The variable angle which this drifting station makes with 
the x-axis is denoted by y. The angle y can take on either positive or negative val- 
ues. As the drifting station moves around the circle, for values of y equal to 
ckj, a n-l' ^ coincides with the 1st, 2d, . . ., (n-l)th stations, respectively. 

Each time such a coincidence occurs the number of stations is effectively reduced by 
one. For this reason, only the situation of n > 2 stations is considered in this part of 
the analysis, as a minimum of two stations must be present for triangulation. 

Substituting equations (A6), (A7), (A8), (A9), and (A 10), which are derived in the 
appendix, into equation (49) leads to 


n-1 

^ ^ip • 2T r j = 5 X 2 cos 2 - 1 + cos 2 yj + 6 y 2 cos 2 0^| + sin 2 yj 
P=0 

+ n6 z ^ sin 2 0 - 26y6 z sin 9 cos 9 sin y 
+ 26^6,, sin 9 cos 9( 1 - cos y) 

+ 6 x 6y cos 2 0 sin(2y) (n > 2) 


(85) 


Substituting equations (48) 
n=l 

for I («p) 2 for n > 2 

P=0 


(50), (51), and (85) into equation (43) results in an expression 
as a function of the elevation angle 9 , the line-of-sight errors 
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d9 and d<j> , the angle a p , and the x-, y-, and z- components of the position error 
P P 


vector 6 r . That is, 


^ (®P f = n ( 6 x 2 + fi y 2 + 6 z 2 ) ' 6 x 2 cos20 (| - 1 + cos2 >) 

p=0 N 

- 6 y 2 cos 2 ^| + sin 2 yj - n6 z 2 sin 2 0 + 2c y 5 z sin 9 cos 9 sin y 

- 26 X 6 Z sin 9 cos 0(1 - cos y) - 6 x 5 y cos 2 0 sin 2 y 

n-1 

+ a 2 ^ |7d0p sin 9 cos 0 !p - d<£ p sin a p J 2 

p=0 ^ 

+ ^d0 p sin 0 sin a p + d^ p cos a p j 2 + ^d0 p cos 0^j 


n-1 


- 2ae x ^ ('dap sin 0 cos a p - d 0 p sin a p \ 
p=o 1 


n-1 


- 2a 6 y ^ ^d0 p sin 0 sin a p + d0 p cos a p j 

p=0 

n-1 


- 2a6 z cos 0 ^ d0 p 

p=0 


(n > 2) 


( 86 ) 


n-1 


In this form, X( 6" p ^ can be minimized with, respect to the x-, y-, and z- components 

P=0 


of the vector 6 r . 


Minimizing equation (86) with respect to 6 X , 6y, and 6 Z , furnishes the following 

three equations, respectively: 

j~n - cos 2 0^| - 1 + cos 2 yJ 6 X + (-cos 2 0 sin y cos y^5 y 

n-1 

+ |sin 0 cos 0(cos y - 1)]6 Z = a Y ^d0 p sin 0 cos a p - d$ p sin oipj (87) 
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(-cos 2 # sin y cos y)6 x + jn - cos 2 #^ + sin 2 yj 6y 

n-1 

+ (sin 8 cos 8 sin y)6 z = a ^ # p sin 8 sin ofp + cos a^j (88) 

P-0 v 


[sin 9 cos 0(cos y - l)]^x + ( s * n ^ coS ® sin y)^y + ( n cos 2 #)6 z = a cos 8 d#p (89) 


Equations (87), (88), and (89) can be rewritten in the following forms, respectively: 


A5 X + B6 y + C6 Z = £ 
B <5 X + D 5y + E § z = Tf 


C6 X + E6 y + F6 Z * $ 


A = n - cos 2 0f — - 1 + cos 2 * 


B = -cos 2 # sin y cos y 


C = sin 0 cos #(cos y - 1) 
D = n - cos 2 #^ + sin 2 ^ 


E = sin 8 cos Q sin y 


F = n cos 2 # 


n-i 

| = a ^ (d0 p sin g cos Ofp - d^p sin ofpj 


41 - 1 

r) = a ^ /d#p sin 8 sin afp + d^p cos ofp\ 

p-0 ' 


( 100 ) 
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and 




= a cos 6 


n-1 

£ % 

p=0 


( 101 ) 


The probLem is now reduced to solving equations (90), (91), and (92) simultaneously 
for 6 X , 6 yJ and 5 Z - Hence, 




B C 
D E 
E F 



- E 2 ) + ?](CE - BF) + £(BE - 



( 102 ) 


A | C 




“y 




= ^r[j(CE - BF) + 7 ]( AF - C 2 ) + £(BC - AE)j 


(103) 


A 



C 


B € 

D T] 

E £ 


= J^Q(BE - CD) + tj(BC - AE) + ?(AD - B 2 )] 


(104) 


where d is the symmetrical determinant 


d- 


ABC 
B D E 

C E F 


(105) 
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Equations (102), (103), and (104) can be rewritten in the following forms, 
respectively: 

6 X = G£ + Htj + I? 

6 y = H£ + 3r} + K£ 

6 Z = I| + Kyj + L? 

where 


(106) 

(107) 

(108) 


G = -^-(DF - E 2 ) 

J&r 

H = -~-(CE - BF) 
& 

I = -^y (BE - CD) 

J = -3j(af - C 2 ) 

K = -i-(BC - AE) 

L = -4(AD - B 2 ) 
jo- 


(109) 

( 110 ) 
( 111 ) 
( 112 ) 

(113) 

(114) 


The coefficients G, H, I, J, K, and L are functions only of n, y, and 9 . 

From equations (106), (107), and (108), respectively, it is recognized that the quan- 
tities 6 X 2 , 6 y 2 , and 6 Z 2 are 

6 X 2 = G 2 ! 2 + hV + I 2 ? 2 + 2GH Jrj + 2Gl|? + 2Hhrf (115) 

5 y 2 = H 2 ! 2 + J 2 ?? 2 + K 2 ? 2 + 2HJ| rj + 2HK?? + 2JKr?? (116) 

6 Z 2 = I 2 | 2 + K V + L 2 ( 2 + 2IK|77 + 2IL|? + 2KL^ (117) 

The quantities | 2 , r) 2 , ij 2 , Jij, f£, and rjif - which are functions, not only of B and 
a, but also of a p , d0 p , and d^> p - need to be evaluated- From equation (99) it can be 
established that the quantity |2 is 
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n-1 n-1 — 

|2 _ a 2 ^ ^ ^d0 p sin 9 cos a p - d<£ p sin a p ^d0 q sin 9 cos a q - d# q sin a q ' 

p=0 q=0 


Using equations (56), (57), and (58) yields 


/ n-1 n-1 \ 

£2 = a2a2jsin20 ^ cos 2o; p + sin^o'pl 

\ P=° p-o / 


Using equations (A 9) and (A 10) renders 


|2 = a 2 o 2 sin 2 - 1 + cos2 y j + + sin 2 y^J 


Using equation (67) results in 


^2 - h 2 a 2 - x + cos2yj + csc 2 0f§ + sin 2 -yj 

(118) 

Similarly, 


^2 = ^2^2^ + sin2yj + esc 2 0^1 - 1 + cos 2 yj 

(119) 

^2 _ h 2 <j2(n cot 2 e) 

(120) 

Jrj = h^o 2 ^-|- cot 2 0 sin 2yj 

(121) 

= h 2 o 2 [cot a(cos y - 1)] 

(122) 

7 }% = h 2 o 2 cot 9 Sin y 

(123) 


If the expressions (118), (119), (120), (121), (122), and (123) are substituted into 
equations (115), (116), and (117) and then these three resulting expressions are substi- 
tuted into equation (65), the rms position error e in the triangulation solution of a point 
in space for the circular configuration of n > 2 observation stations with equal elevation 
angles, when one of the stations is permitted to drift around the circle, can be computed. 
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The expression for e/ha for n = 3 was programed for a high-speed computer to get 
numerical values. 

The quantity e/ha was evaluated for values of y - 0° to 360°, at 5° intervals, for 
each of the values of 6 = 5° to 85°, at 5° intervals. For each particular value of 6 
examined, the quantity e/ha was the smallest when y ~ 0° (and 360°). Table HI is 
a table of values of e/ha as a function of both 9 and y. For convenience, the values 
listed in the table are only those for 9 and y at 20° intervals. It is observed from 
table HI that for each particular value of 6 , the smallest value of e/ha occurs when 
y = 0° (and 360°). 

Figure 5 is a plot of e/ha as a function of y for 9 = 5° to 85°. Again, for con- 
venience, 20° intervals in 9 were used. It is again seen from figure 5 that for each 
particular value of 6 the quantity e/ha is a minimum when y = 0° (and 360°), which 
represents the symmetrical arrangement of the circular configuration of observation sta- 
tions with equal elevation angles. 

It can be seen from table III and figure 5 that when y = 180° (i.e., when the drifting 
observation station is again equidistant from the two fixed stations), then e/ha is either 
a relative minimum or a maximum, depending on the value of the elevation angle 9. It 
is a relative minimum for 6 equal to 5° and 25°; and it is a maximum for 9 equal to 
45°, 65°, and 85°. For 9 equal to 5° and 25°, the quantity e/ha is a maximum when 
y is equal to 90° and 270°. 

Figure 6 is a plot of e/ha as a function of the elevation angle 9 for y = 0°, 60°, 
120°, and 180°. The quantity e/ha was examined for y only through 180°, as e/ha 
is symmetrical with respect to y - 180°, which can be seen from table III. It is observed 
from figure 6 that for each particular value of y the quantity e/ha is a minimum for 
6 = 55.300°. It is remembered that for the symmetrical circular configuration of n > 2 
observation stations with equal elevation angles the optimum elevation angle, which min- 
imizes the rms position error in the triangulation solution of a point in space, is 55.300°. 

Therefore, the symmetrical arrangement for n = 3 observation stations with equal 
elevation angles in a circular configuration giveg the minimum rms position error in the 
triangulation solution of a point in space. It may be plausibly assumed, then, that the 
symmetrical arrangement is also the best arrangement for n > 3 stations with equal 
elevation angles in a circular configuration. 

CONCLUSIONS 

The purpose of this paper was to derive, to first order, analytical expressions for 
the rms position error in the triangulation solution of a point object in space for several 
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ideal observation- station configurations to provide insights into the nature of the depend- 
ence of the rms position error on certain of the experimental parameters involved. 

For two arbitrarily located observation stations, the rms position error in the tri- 
angulation solution of a point in space is directly proportional to the rms line -of -sight 
error. The minimum rms position error occurs when the two elevation angles are both 
equal to 56.098°. For the configuration of the two arbitrarily located stations, two choices 
of the most probable point in space - the midpoint of the shortest line between the two 
lines of sight, and the point, on the aforementioned line, which subtends equal residual 
angles at the two stations - were used. The latter choice of the most probable point 
yields the smaller value of the rms position error, except in the special situation of 
equal elevation angles and then the two choices yield the same value. For the first choice 
of the most probable point, the optimum relative weighting factor, which minimizes the 
rms position error, for the data from the first station when the data from one of the sta- 
tions is degraded relative to that from the other, is directly proportional to the square of 
the ratio of the rms line-of-sight errors from the second and first stations, respectively. 

For the symmetrical circular configuration where the number n of observation 
stations with equal elevation angles is two or more, the rms position error in the trian- 
gulation solution of a point in space is directly proportional to the rms line-of-sight 
error. For n > 2 the rms position error is inversely proportional to the square root 
of the total number of stations of interest. For n = 2 the optimum elevation angle, 
which minimizes the rms position error, is 56.098°. For elevation angles of 56.098° ± 5°, 
the difference in the rms position error for n » 2 is 1.5 percent. For n > 2 the 
optimum elevation angle, which minimizes the rms position error, is 55.300° regard- 
less of how many (greater than two) stations are present. For elevation angles of 
55.300° ± 5°, the percent difference in the rms position error for n - 3 is 2.5 percent, 
while for n = 8 it is 0.9 percent. The value of the optimum elevation angle for n > 2 
is within a fraction of a degree of the angle (54.667°) which the sloping edge of a tetra- 
hedron makes with its base. Therefore, the optimum configuration of three stations for 
triangulating on a point in space is the configuration of the three stations plus the point 
which closely approximates a tetrahedron. 

When one of the stations in the circular configuration of n = 3 observation stations 
with equal elevation angles is permitted to drift around the circle from its position of 
symmetry, the rms position error in the triangulation solution of a point in space is a 
minimum when the variable angle of the drifting station equals zero degrees, which repre- 
sents the symmetrical arrangement. For each value of the variable angle of the drifting 
station, the rms position error is a minimum when the elevation angle equals 55.300°. 

Based on the assumptions used in this study, a necessary and sufficient condition 
for a minimum rms position error in the triangulation solution of a point in space for two 
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arbitrarily located observation stations is that the two elevation angles be equal, and equal 
to 56.098°. A necessary and sufficient condition for a minimum rms position error in 
the solution of a point for three stations with equal elevation angles in a circular con- 
figuration is that the three stations be symmetrically arranged around the circle. It 
may be plausibly assumed, then, that the symmetrical arrangement of the circular con- 
figuration of n > 3 stations -with -equal elevation angles is the arrangement of this con- 
figuration which yields the minimum rms position error. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., March 14, 1974. 
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APPENDIX 


DERIVATION OF CERTAIN TRIGONOMETRIC SERIES USED FOR 
CIRCULAR CONFIGURATION OF n g 2 STATIONS 
WITH EQUAL ELEVATION ANGLES 


For the symmetrical circular configuration of n g 2 observation stations with 
equal elevation angles, the angle between the x-axis and the pth station is a p = 2irp/n, 


as shown in figure 2. Now, 


n-1 

n-1 

^ cos a p = Re 

i exp (‘ — ) 

p~o 

A \ t 

P“0 


= Re 


— 

exp(i2v) 

- 1 

exp 


- 1 


Therefore, 


n-1 

i 

P=0 


£ COS 0!p = 0 


Similarly, 


n-1 

i 

P=0 


sin Q!p » 0 


Also, 


n-1 

sin 2a p - Im 

p=0 


n-1 

i 

P-0 


exp i 


4?rp 


r 


ss Im 


exp(i47r - 1) 


ex p( 1 ¥ - l ) 


For n = 2, 


(n g 2) (Al) 


(n g 2) (A2) 


0 (n > 2) 


n-1 


^ sin 2a p = ^ sin(2^p) = 0 

p=0 p=0 
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APPENDIX - Continued 


Therefore, 


n-1 

p=0 


sin 2ap = 0 


(n 2) (A3) 


In addition, 


n-I 

P=0 


cos 2a p = Re| 


n-1 


£ exp(i 4 * p 


P=0 


= Re 


exp(i47r - 1) 


exp(i ^ - 1 




= 0 (n > 2) 


For n = 2, 


n-1 


^ cos 2a p = cos(2itp) = 2 

p=0 p=0 


Hence, 


n-1 

/ cos 2a D = 

p=0 12 


(n >2) 
(n = 2) 


Since 


n-1 


n-1 

) cos^a« 

_ 1 

y 

L j P 

2 

L 

p=0 


p=0 


1 + cos 2a^J 


then 


n-1 ( n 

,2„ J* 


I 

p=o 


cos*a. 


(n > 2) 
(n = 2) 


(A4) 
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APPENDIX - Continued 


Since 


then 


n-1 


n-1 


I sin2 «p = II ( x - cos 2 “p) 

p=0 p=0 


n-1 

^ sln2 “ P = 

p=0 



(n>2) 
(n = 2) 


(A 5) 


If, for n > 2 observation stations, the Oth station j^i.e., the one previously located 
on the x-axis and for which otp = Q!q = 0^ is permitted to drift around the circle, making 
the variable angle y with the x-axis, then equations (Al), (A2), (A3), (A4), and (A5) must 
be modified- Now, 


n-1 

p-v 


n-1 


COS Q!p = COS Q!q + 


cos a, 


p~l 


Therefore, 


n- j 

l 


P -0 


COS Q!p 


= cos y - 


1 


Also, 


n-1 n=l 

^ sin otp = sin «g + ^ s * n a p 
p=0 p=l 


Therefore, 



sin afp = sin y 


In addition, 


n- j 

I 


P=0 


sin 2afp = sin 2ai g + 


n- j 

I 


p=l 


sin 2a p 


(n > 2) (A6) 


(n > 2) (A7) 
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APPENDIX - Concluded 


Therefore, 


n-1 

sin 2ofp = sin 2 y 

p=0 


Furthermore, 

n-1 


n-1 


^ COS^Qip = COS2 q!q + ^ COS^CKj 


Therefore, 


P-0 


n-1 


p=l 


Finally, 


Therefore, 


^ COS^QJp ■ COS^y + | - 1 
p=0 


n-1 n-1 

sin2o!p = sin^aiQ + ^ sin2a>p 

p=0 p=l 


n-1 

^ sin2ap = sin^y + H 
P=0 


(n > 2) (A8) 


(n > 2) (A 9) 


(n > 2) (A 10) 
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TABLE I.- THE QUANTITY e/ha IN SOLUTION OF POINT FOR TWO ARBITRARILY 
LOCATED STATIONS, WITH MOST PROBABLE POINT TAKEN AS MIDPOINT OF 
SHORTEST LINE BETWEEN LINES OF SIGHT, AS A FUNCTION OF THE 
TWO ELEVATION ANGLES IN INTERVALS OF 10° 


0 

0 2 , deg 

— , — 

deg 

5 

15 

25 

35 

45 

55 

65 

75 

85 

5 

93.79 

35.91 

24.15 

18.96 

16.16 

14.52 

13.55 

13.04 

12.88 

15 

35.91 

11.26 

7.40 

5.93 

5.18 

4.76 

4.55 

4.47 

4.52 

25 

24.15 

7.40 

4.68 

3.70 

3.24 

3.01 

2.92 

2.92 

3.02 

35 

18.96 

5.93 

3.70 

2.90 

2.54 

2.38 

2.34 

2.38 

2.53 

45 

16.16 

5.18 

3.24 

2.54 

2.24 

2.11 

2.11 

2.21 

2.42 

55 

14.52 

4.76 

3.01 

2.38 

2.11 

2.03 

2.07 

2.24 

2.58 

65 

13.55 

4.55 

2.92 

2.34 

2.11 

2.07 

2.18 

2.47 

3.08 

75 

13.04 

4.47 

2.92 

2.38 

2.21 

2.24 

2.47 

3.02 

4.28 

85 

12.88 

4.52 

3.02 

2.53 

2.42 

2.58 

3.08 

4.28 

8.21 
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TABLE H.- THE QUANTITY e/ho IN SOLUTION OF POINT FOR TWO ARBITRARILY 
LOCATED STATIONS, WITH MOST PROBABLE POINT TAKEN AS MIDPOINT OF 
SHORTEST LINE BETWEEN LINES OF SIGHT, AS A FUNCTION OF THE 
TWO ELEVATION ANGLES IN INCREMENTS OF 0.001° 


*v 

deg 

e 2 , deg 

56.096 

56.097 

56.098 

56.099 

56.100 

56.096 

2.0278919610 

2.0278919578 

2.0278919558 

2.0278919550 

2.0278919554 

56.097 

2.0278919578 

2.0278919550 

2.0278919535 

2.0278919532 

2.0278919541 

56.098 

2.0278919558 

2.0278919535 

2.0278919525 

2.0278919527 

2.0278919540 

56.099 

2.0278919550 

2.0278919532 

2.0278919527 

2.0278919533 

2.0278919552 

56.100 

2.0278919554 

2.0278919541 



2.0278919540 

2.0278919552 

2.0278919575 
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TABLE IQ.- THE QUANTITY e/ha IN SOLUTION OF POINT FOR CIRCULAR 
CONFIGURATION OF n = 3 STATIONS WITH EQUAL ELEVATION ANGLES, 
WHEN ONE STATION IS PERMITTED TO DRIFT AROUND 
THE CIRCLE, AS A FUNCTION OF THE ELEVATION 
ANGLE AND THE VARIABLE ANGLE 


y, 

0, deg 

. 

deg 

5 

25 

45 

65 

85 

0 

14.77913 

2.93361 

1.76383 

1.77858 

6.70017 

20 

15.08990 

2.96540 

1.77248 

1.78803 

6.74501 

40 

16.01071 

3.05499 

1.79750 

1.81630 

6.87938 

60 

17.38405 

3.17939 

1.83586 

1.86287 

7.10196 

80 

18.56983 

3.28970 

1.88192 

1.92612 

7.40728 

100 

18.58380 

3.32772 

1.92795 

2.00217 

7.78023 

120 

17.41489 

3.27762 

1.96638 

2.08354 

8.18797 

140 

16.04393 

3.18378 

1.99285 

2.15846 

8.57282 

160 

15.12098 

3.10440 

2.00725 

2.21226 

8.85561 

180 

14.80908 

3.07453 

2.01166 

2.23197 

8.96061 

200 

15.12098 

3.10440 

2.00725 

2.21226 

8.85561 

220 

16.04393 

3.18378 

1.99285 

2.15846 

8.57282 

240 

17.41489 

3.27762 

1.96638 

2.08354 

8.18797 

260 

18.58380 

3.32772 

1.92795 

2.00217 

7.78023 

280 

18.56983 

3.28970 

1.88192 

1.92612 

7.40728 

300 

17.38405 

3.17939 

1.83586 

1.86287 

7.10196 

320 

16.01071 

3.05499 

1.79750 

1.81630 

6.87938 

340 

15.08990 

2.96540 

1.77248 

1.78803 

6.74501 

360 

14.77913 

2.93361 
-- — 

1.76383 

1.77858 

6.70017 
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E/ho 



0,deg 

Figure 4.- The quantity e/ha in solution of point for symmetrical circular 
configuration of n = 2 stations with equal elevation angles as a function 
of the elevation angle . 
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Figure 5.- The quantity e/ha in solution of point for circular configuration 
of n = 3 stations with equal elevation angles, when one station is per- 
mitted to drift around the circle, as a function of the variable angle for 
different values of the elevation angle. 
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Figure 6.- The quantity e/hcr in solution of point for circular configuration 
of n = 3 stations with equal elevation angles, when one station is per- 
mitted to drift around the circle, as a function of the elevation angle for 
different values of the variable angle. 
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